Strangeness equilibration in heavy ion collisions 
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Using a relativistic transport model for heavy ion collisions at energies that are below the thresh- 
old for kaon and antikaon production in nucleon-nucleon collisions, we study how their abundances 
approach the canonical equilibrium during the collisions. We find that kaons are far from chemical 
equilibrium at the initial and high density stage, and they approach equilibrium only during the ex- 
pansion stage of the collisions when their production rate is small and becomes comparable to their 
annihilation rate. In contrast, antikaons approach chemical equilibrium much earlier but eventually 
fall out of equilibrium again as a result of their large annihilation cross sections in nuclear matter. 
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Recent analyses have shown that most hadrons mea- 
sured in heavy ion collisions can be described by statis- 
tical models based on the grand canonical ensemble for 
abundant particles ^ and the canonical ensemble for rare 
particles [0. On the other hand, studies based on trans- 
port model indicate that the chemical equilibration time 
for kaons and antikaons in the hot dense matter, that 
is expected to be formed in heavy ion collisions at the 
' AGS and SPS energies, is an order of magnitude longer 
than the heavy ion collision time 1^. It has thus been 
suggested Q that the kaon equilibration time can be sig- 
nificantly shortened if the kaon mass is reduced in dense 
matter as a result of the large attractive scalar interac- 
tion and diminishing repulsive vector interaction due to 
the assumption of vector decoupling. The disappearance 
of the repulsive vector interaction is, however, not con- 
sistent with the strong repulsive vector potential that is 
' needed to understand the vanishing kaon flow observed 
in heavy ion collisions at GSI energies [HJ|] and the large 
kaon antiflow observed in heavy ion collisions at AGS 
energies 0,1). Also, a recent transport model study has 
shown that the measured kaon yield in heavy ion colli- 
sions at GSI 1^ can be explained with kaons interacting 
with both the scalar and vector potentials . 

To gain insight into this problem, a kinetic theory 
for the time evolution of particle production pl] | can 
be applied. It is found that the equilibrium time, Tq, 
for rare particles carrying U{1) charge and described by 
the canonical ensemble due to U{1) charge conservation, 
is much shorter than what is expected from the grand 
canonical ensemble, i.e., 
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where is the average multiplicity of rare particles 

per event if they were described by the grand canonical 
ensemble. Also, it is shown that 
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so the equilibrium multiplicity in the canonical ensem- 
ble is much lower than that given by the grand canonical 
ensemble. 

The above idea based on the kinetic theory can be 
quantitatively studied using the transport models, which 



have been very successful in understanding many aspects 
of heavy ion collisions |13|. In this paper, we shall use 
the transport model to study kaon and antikaon produc- 
tion in heavy ion collisions at energies that are below the 
thresholds for their production in nucleon-nucleon inter- 
actions, and to understand how kaons and antikaons ap- 
proach chemical equilibrium in these collisions. In previ- 
ous studies of subthreshold kaon production using trans- 
port models ||l^,|l^, only the production of kaons has 
been included while their annihilation via the inverse re- 
actions has been neglected. The neglect of kaon anni- 
hilation would not be a bad approximation if kaons are 
far from chemical equilibrium. However, to understand 
kaon chemical equilibration in transport models, we need 
to include kaon annihilation to see how the kaon produc- 
tion and annihilation rates become comparable during 
the finite heavy ion collision time. 

We follow essentially the same model used in Ref. , 
which is based on the relativistic transport model RVUU 
pTf . In this model, the nuclear potential is taken from 
the nonlinear Walecka model, so it has both an attractive 
scalar and a repulsive vector part. The attractive scalar 
potential allows one to treat consistently the change of 
nucleon mass in nuclear matter. In dense matter, the 
nucleon mass is reduced and the energy is in the scalar 
field. As the system expands, the nucleon regain its mass 
from the scalar field energy. The model in Ref. ||l6| in- 
cludes kaon production from both baryon-baryon and 
meson-baryon interactions. For kaon production from 
baryon-baryon interactions, the cross sections are taken 
from the predictions of the meson-exchange model intro- 
duced in Ref. |l^. For kaon production cross sections 
from meson-baryon interactions, they are obtained from 
the resonance model of Ref. |1^. The produced kaons 
together with their partners, mainly hyperons, not only 
undergo elastic scatterings with baryons but are also af- 
fected by mean-field potentials. For kaon and hyperon 
scattering cross sections, we take the empirical values as 
in Ref. [ |l6| . The kaon potential is taken from the chiral 
Lagrangian including both scalar and vector interactions 
pO| with their strengths determined from experimental 
observables such as the kaon yield and collective fiow in 
heavy ion collisions U&. 
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Since the kaon production probability is much smaller 
than one in heavy ion collisions at subthreshold ener- 
gies, it is usually treated by the perturbative method, 
i.e., the effect of kaon production on the dynamics of 
heavy ion collisions is neglected. Specifically, a kaon 
is produced in a baryon-baryon or a meson-baryon in- 
teraction that is above the production threshold. The 
produced kaon then carries a probability given by the 
ratio of the kaon production cross section to the total 
baryon-baryon or meson-baryon cross section. Further- 
more, to treat the rescattering of kaons and hypcrons 
with other particles, one allows more than one pair of 
kaon and hyperon to be produced in such a collision. 
When these kaons and hyperons are scattered by other 
particles, their momenta are changed according to the 
differential cross section, which is usually taken to be 
isotropic, while changes of the momenta of the other par- 
ticles are neglected. The final kaon yield and spectrum 
are then obtained by adding the probabilities of individ- 
ual kaons and dividing by the number of kaons produced 
in each baryon-baryon or meson-baryon collision. 

To study the chemical equilibration of kaons, we im- 
prove the model of Ref . [|6j by including also kaon anni- 
hilation by hyperon, using the cross sections determined 
from the meson-baryon interactions via the detailed bal- 
ance relations. Because of strangeness conservation, a 
kaon is produced together with a hyperon. As the kaon 
production probability is much less than one in heavy 
ion collisions at subthreshold energies, there is only one 
hyperon in an event in which a kaon is produced. Kaon 
annihilation can thus occur only when there is a colli- 
sion between the same pair of kaon and hyperon that is 
produced in the baryon-baryon or meson-baryon inter- 
action. As shown in Ref. [pl| via the kinetic equation, 
the annihilation between such a pair of particles that are 
produced simultaneously as a result of the U{1) charge 
conservation would lead to an equilibration described by 
the canonical ensemble. 

To illustrate the physics of kaon chemical equilibra- 
tion in heavy ion collisions, we consider Ni-|-Ni collisions 
at lA GeV and impact parameter 6 = fm, which is 
below the threshold for both kaon and antikaon produc- 
tion in nucleon-nucleon interactions. Similar results are 
obtained for collision energies at 2A GeV, which is above 
the kaon production threshold but below the antikaon 
production threshold in nucleon-nucleon interactions. In 
the top panel of Fig. ^ the time evolution of the kaon 
abundance for the scenarios with (solid curve) and with- 
out (dotted curve) kaon annihilation are given. It is seen 
that including kaon annihilation by hyperon reduces the 
final kaon yield by only about 10%. As shown in the 
lower window of Fig. |l|, the kaon production rate (solid 
curve) is appreciable only when the nuclear density (thick 
solid curve) is high. The effect due to kaon annihilation 
is better illustrated by the kaon annihilation rate (dashed 
curve) shown in the bottom panel of Fig. |l|. We see that 
the annihilation rate is negligible during the high den- 
sity stage when most kaons are produced. This result 



thus justifies the neglect of kaon annihilation in previous 
studies. Although the kaon annihilation rate is small, 
it becomes comparable to the kaon production rate at 
about 17 fm/c, indicating that the kaon yield eventu- 
ally approaches chemical equilibrium. We note that the 
baryon density at which kaons reach chemical equilib- 
rium is about 1.2 po, where po ~ 0.16 fm~^ is the normal 
nuclear matter density. 
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FIG. 1. Top panel: Time evolution of kaon abundance 
in Ni-|-Ni collisions at lA GeV and impact parameter 6 = 0. 
Dotted curve is the result without kaon absorption while solid 
and dashed curves are results with kaon absorption from us- 
ing default KB and YB elastic scattering cross sections and 
twice the cross sections, respectively. Bottom panel: Time 
evolution of kaon production rate (solid curve) and absorp- 
tion rate (dashed curve) as well as the central baryon density 
(thick solid curve). 



To better understand kaon production in heavy ion 
collisions, it is useful to know if the system is close to 
thermal equilibrium when kaons are produced. For this 
purpose, we have evaluated the average squared momen- 
tum of all particles both in the beam direction, (pi), and 
in the transverse direction, (p^). The ratio {pl)/{Px) can 
then be used to characterize the degree of thermal equi- 
librium, with a value of one corresponding to complete 
thermal equilibrium. In Fig. H, we show the time evolu- 
tion of this ratio for nucleons (solid curve), deltas (dotted 
curve), and pions (dashed curve) in the central volume 
of 8 fm'^. It is seen that at the time when most kaons 
are produced, both pions and deltas are close to ther- 
mal equilibrium while nucleons are not. Since kaons are 
dominantly produced from pion and nucleon interactions 



2 



(about 70%) , our results thus show that they are mostly 
produced from a dense but not thermally equilibrated 
hadronic matter. 
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FIG. 2. Time evolution of the ratio {pi) / {pi) for nucleons 
(solid curve), deltas (dotted curve), and pious (dashed curve). 
The central baryon density as a function of time is given by 
the thick solid curve. 

According to the kinetic theory , the chemical equi- 
libration time in the canonical formalism is given by 
T§ = V/L = NK/idNL/dt), where V is the volume 
of the region where kaon annihilation occurs, L is the 
momentum averaged cross section for kaon annihilation, 
and Nk and dN^/dt are the kaon number and its ab- 
sorption rate, respectively. At time t = 12 fm/c when 
the kaon absorption rate is largest, we have from Fig. |^, 
Nk = 4.2 X 10-3 and dNL/dt = 7.0 x 10-^ which give 
a kaon chemical equilibration time of 60 fm/c if the sys- 
tem is prevented from expanding from t = 12 fm/c. On 
the other hand, the thermal average of the kaon annihila- 
tion cross section (Tky^-kN is about 0.25 fm^ at temper- 
ature T = 75 MeV, and it changes by less than 20% for 
50 < r < 100 MeV due to the exothermic nature of the 
annihilation process. The above value for the kaon chem- 
ical equilibration time implies that the effective volume 
in which kaon annihilation occurs is about 15 fm'^. Since 
the chemical equilibration time for kaons is much longer 
than the heavy ion collision time, we would normally ex- 
pect them not to reach chemical equilibrium during the 
collisions. However, since the kaon production rate de- 
creases strongly as the temperature decreases due to the 
large threshold of the production process, its value can 
thus become comparable to the annihilation rate when 
the system expands and cools. When this happens for 
later stage of heavy ion collisions as shown in the trans- 
port model results of Fig. |l|, kaons can then reach chem- 
ical equilibrium. We thus note that, since most kaons 
are produced in the dense but non-equilibrium stage, the 
chemical equilibration time could be overestimated when 
determined in the normal way. Furthermore, the chemi- 
cal freeze-out temperature extracted from the final kaon 



abundance in thermal model analysis does not provide 
sufficient information on the dynamics of kaon produc- 
tion. 

The effect due to kaon annihilation depends on the 
magnitude of kaon and hyperon elastic scattering cross 
sections with other particles. If there are no such scat- 
terings, e.g, if these cross sections are set to zero, then 
the produced kaon and hyperon would simply move away 
from each other without further interactions, leading to a 
result similar to that without kaon annihilation. On the 
other hand, larger kaon and hyperon scattering cross sec- 
tions with other particles would force them into a smaller 
region, thus increasing the kaon annihilation rate. This 
is demonstrated in the top panel of Fig. |l| by the dashed 
curve, which is obtained by taking the kaon and hyperon 
scattering cross sections with other baryons to be twice 
the default values. It is seen that these larger cross sec- 
tions indeed further reduce the kaon yield. 
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FIG. 3. Same as Fig. |l| for antikaons. 

We also study antikaon production within the relativis- 
tic transport model, where antikaons are produced not 
only from baryon-baryon and meson-baryon interactions 
but also from meson-hyperon interactions. The cross sec- 
tions for antikaon production from both baryon-baryon 
and meson-baryon interactions are taken from an anal- 
ysis based on the meson-exchange model |2|]. For an- 
tikaon production from the meson-hyperon interactions, 
the cross sections are obtained from the empirical cross 
sections for K~ absorption by nucleon to form a pion 
and a hyperon |^2|. As first pointed out in Ref. |23|, an- 
tikaon production in heavy ion collisions at subthreshold 
energies is mainly due to the meson-hyperon interactions. 
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For antikaon annihilation, it is dominated by the reaction 
KN — > -kY , and this has aheady been included in previ- 
ous transport model studies [ p^|j2^ , |25| ] . We note that to 
account for the observed enhancement of antikaon pro- 
duction in heavy ion collisions at subthreshold energies, 
it has been shown in Refs. |l^j2j,^ that a dropping of 
antikaon mass due to medium effects is needed. 

The results for antikaon production in Ni-fNi collisions 
at lA GeV and impact parameter 6 = fm are shown in 
Fig. ^. In the top panel, the time evolution of antikaon 
abundance is shown for the cases with (solid curve for de- 
fault kaon and hyperon elastic scattering cross sections 
and dashed curve for twice the default cross sections) and 
without (dotted curve) kaon annihilation. The reduction 
of antikaon yield when kaon annihilation is allowed is 
due to the reduction in the production probability of hy- 
perons, which contribute most to antikaon production. 
The time evolution of the antikaon production and an- 
nihilation rates are shown in the bottom panel of Fig. ^ 
by the solid and dashed curves, respectively. Similar to 
kaon production, we find that antikaons are mostly pro- 
duced in the high density stage of heavy ion collisions. 
Because of the large baryon density, their annihilation 
rate through the reaction KN ttY is larger than the 
kaon annihilation rate. As a result, antikaons approach 
chemical equilibrium even in the earlier stage of heavy 
ion collisions. However, they eventually fall out of equi- 
librium as shown in the lower panel of Fig. ^ 

To summarize, we have improved the perturbative 
treatment of strange particle production in the rela- 
tivistic transport model by including the annihilation of 
kaon and hyperon that are produced as a pair due to 
strangeness conservation. For heavy ion collisions at en- 
ergies below the threshold for strange particle produc- 
tion in nucleon-nucleon interactions, we find that both 
kaons and antikaons are largely produced during the high 
density stage of the collisions when the system has not 
reached thermal equilibrium. Because of their large anni- 
hilation cross sections in dense nuclear matter, antikaons 
are near chemical equilibrium much earlier. For kaons, 
their abundance at high density stage is far from equilib- 
rium, and it only becomes close to the equilibrium value 
during the expansion stage of heavy ion collisions when 
the production rate is small and comparable to the anni- 
hilation rate. The small kaon annihilation rate in heavy 
ion collisions thus justifies the neglect of kaon annihila- 
tion in previous transport model studies of subthreshold 
kaon production. 
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